EP 924 910 



Page 1 of 26 



(19) 




EuropSisches Patentamt 
European Patent Office 
Office europden des brevets 




(11) 



EP0 924 910 A1 



(12) 



EUROPEAN PATENT APPLICATION 



(43) Date o1 publication: 

23.06.1999 Bulletin 1999/25 

(21) Application number: 98310269.0 

(22) Dated filing: 15.12.1998 



(51) Int Cl.^: H04L 27/26, H04L 1/00. 

H04L 25/49. H04L5/06 



(84) 


Designated Contracting States: 


• Davis, James Andrew 




AT BE CH CY DE DK ES R FR GB GR IE IT LI LU 


Richmond, VA 23233 (US) 




MC NL PT SE 


• Peterson, Kenneth Graham 




Designated Extension States: . 


, Cotham, Bristol BS6 5TP (GB) 




ALLTLVMKROSI 


(74) Representative: 






(30) 


Priority: 1 7.1 2.1 997 EP 9731 0252 


Coker, David Graeme et al 




Hewlett-Packard Limited 


(71) 


Applicant: 


intellectual Property Section 




Hewlett-Paclcard Company 


Building 2 




Palo Alto, California 94304 (US) 


RIton Road 






Stoke Gifford, Bristol BS34 8QZ (GB) 


(72) 


Inventors: 




Jedwab, Jonathan 

Oalcf ieid Road, Clifton, Bristol BS8 2BA (GB) 





(54) Method and apparatus for decoding COFDii/l-data 

(57) In a coded orthogonal frequency division multi- 
plex (COFDM) systenr /7-bit data words are enc(xled as 
2 '"-symbol code words, each symbol having 2^ possible 
values (e.g. y = 3 for octary). An efficient decoder for 
these code words is provided by applying j iterations of 
the fast Hadamard transform, the input vector for the 
second and subsequent iterations being derived from 
the result of the immediately preceding iteration. The 
decoder is also applicable to p^-ary codes in which the 
code words symbols each have pi possible values, 
where p is a prime other tiian 2. 
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Description 

4 

Technical Field 

5 [0001] This invention relates to methods and apparatus for decoding data, and to communications systems incorpo- 
rating such methods and apparatus; in particular, though not exclusively, the invention relates to such methods, appa- 
ratus and systems for decoding data transmitted using multi-carrier frequency division multiplexing, such as coded 
orthogonal frequency division multiplexing (COFDM). 

10 Background Art 

[0002] The increasingly widespread use of powerful computer-lsased equipment and appliances is resulting in rapidly 
increasing requirements for high-speed, high-capacity communications, based on both wireline and wireless tech- 
niques. One way of accommodating this inaease is to design communications systems for operation at higher and 

75 higher symbol rates. However, an increase in symbol rate typically involves a corresponding decrease in inter-symbol 
spacing, and this in turn can lead to greater inter-symbol interference (ISl) and consequential errors in symbols as 
received. In the case of wireline links ISl can arise, for example, from signal reflections at imperfectly matched inter- 
faces between different portions of the communications link. For wireless links a significant source of iSI is multi-path 
propagation involving signal reflections from objects (such as buildings) in the space traversed by the signals; mobile 

20 wireless communications can be especially susceptible to ISl from this source, because movement of a mobile unit 
using the link can cause continuous and complex variation In the nujltiple propagation paths extending to the mobile 
unit. 

[0003] A technique which is being actively studied and developed to Increase communications system capacity and 
counteract ISl is multi-camer (or multitone) frequency division multiplex operation, such as COFDM operation. In a 

25 COFDM system multiple symbols are transmitted simultaneously, each symbol controlling the modulation (e.g. by 
phase-shift Keying) of a respective one of multiple carrier signals having closely-spaced frequencies. By choosing the 
can'ier-frequency spacing to have an appropriate relationship to the rate at whidi symbols are transmitted on each car- 
rier, it is possible to ensure that each individual carrier can be demodulated without interference from the carriers adja- 
cent to it (in which case the carriers are said to be orthogonal to one another). The set of symbols transmitted 

30 simultaneously on respective carriers during a symbol period (a code word) is chosen to encode a group of data sym- 
bols (a data word) and includes redundancy (i.e. the code word has more symtx)ls than the data word); this redundancy 
is used upon reception for detection of errors in the symbols as received, and. if desired, for correction of those errors. 
[0004] COFDM systems counteract ISl in part because the transmission of multiple symbols in parallel enables a 
desired data rate to be attained with a lower symbol rate than would be necessary if a serial system were used to trans- 

35 mit individual symbols In succession. A lower symbol rate Implies a longer duration for each symbol and thus a longer 
period between symbol transitions, which reduces the impact of ISl. The error detection/carectlon capabilities provided 
by the inclusion of redundancy enhance the resistance to data con-uption, by enabling some errors which occur in 
received symbols owing to signal fading and other causes to be deteqted or corrected. 

[0005] However, present COFDM systems encounter difficulties in obtaining the maximum potential benefit for a par- 
40 ticular transmitted signal power, owing to limitations arising from the typical peak-to-mean envelope power ratio 
(PMEPR) of the complex-valued signal from which the actual transmitted signal is derived. Because this transmitted 
signal is effectively the sum of several signals at closely spaced frequencies and with respective phases which change 
at each symbol transition, the overall signal tends to exhibit transient pronounced peaks separated by substantially 
longer intervals of much lower though still varying amplitude. The transmitting equipment must be adjusted to reproduce 
45 the peaks without dipping or other distortion, so for much of the time, in the intervals between peaks, that equipment is 
operating at a signal level well below its maximum capability. As a result the geographical range of the equipment is sig- 
nificantly lower than would be expected for its nominal power rating, or alternatively a much more powerful equipment 
must be provided to cover the desired range and then operated inefficiently 

[0006] It is known that this problem can in principle be alleviated by careful selection of the code words which encode 
so the possible data words (for example in 'Block coding scheme for reduction of peak to mean envelope power ratio of 
mutticarrier transmission schemes' by A.E. Jones, T.A. Wilkinson and S.K. Barton, Electronics Letters, 8tii December 
1994. vol.30, no.25, pp. 2098-2099). However, in selecting code words to reduce tiie PMEPR it is important also to 
ensure tiie code has good error correction properties. Furthermore, implementation of a high-tiiroughput comnmjnica- 
tions system desirably involves the use of long code words and of rapid encoding and decoding procedures; the use of 
55 simple look-up tables, though readily implemented, does not lend Itself to rapid operation using compact (and low- 
power, inexpensive) circuitry, espedalty during decoding, if the number of valid code words to be checked Is large, as 
is typically the case with long code words (e.g. more tiian about 16 symbols). It is therefore preferable for encoding and 
decoding operations to be definable in terms of some analytical procedure which for example performs combinatorial 



2 



£P&'24 910 



Page 3 of 26 



EP 0 924 910 A1 

logic on the data word or code word symbols to corrvert one to the other, or at least restricts the size of any look-up table 
which is used. 

[0007] An addrtionai difficulty is presented by the sheer quantity of possible code words from which a selection may 
be made to define a code for use in a typical system. As the number of symbols in a code word inaeases the number 

5 of possible code words increases much faster, even for binary symbols; for symbols having just three (ternary) or four 
(quaternary) possible values, the number of possible codes available becomes so huge even for quite short code words 
that a simple exhaustive search through ail possible code words is not feasible in a realistic time even just to minimize 
the PMEPR. Attempting to satisfy enror connection and implementational requirements as well merely adds to the diffi- 
culty of such a search, especially as it will not even be known whether there actually exists a code word set which is of 

10 the desired size and which simultaneously satisfies prefen-ed criteria relating to PMEPR, error detection/correction and 
ease of encoding and decoding. This difficutty is further exacerbated by the desirability of having a large number of 
code words available for use in the chosen code, in order to be able to encode a large number of different data words 
■ (i.e. data words containing relatively many symbols) and thus attain an acceptably high data transfer rate. 

[0008] It is known that certain kinds of codes may have some desirable properties. Thus the paper 'Synthesis of 

75 Power Efficient Multltone Signals with Rat Amplitude Spectrum' by B.M. Popovi6, IEEE Transactions on 
Communications, vol.39, no.7, July 1991, pp. 1031 -1033 indicates that any binary or polyphase complementary 
sequence can be used to construct a multitone signal such that the crest factor of the real-valued signal is less than or 
equal to 6 dBV (which corresponds to the PMEPR of the complex-valued signal being less than or equal to 3 dB). Exam- 
ples of these sequences have been described, as in the case of binary Qolay complementary sequence pairs ('Com- 

20 plementary series' by M.J.E. Golay, IRE Transactions on Information Theory, vbl.lT-7, Apr. 1961, pp. 82-87). However. 
Popovi6's paper is directed to situations such as multipath fading measurements and multitone jamming, in which a sin- 
gle code word can be used continuously In generating the required signal. Acoordingty Popovi6 gives no guidance on 
how multiple such sequences could be identified just to attain the 3 dB limit on PMEPR in a practical code for conveying 
information; and no mention whatsoever is made of the other requirements mentioned above for a code suitable for 

25 practical implementation. Indeed Popovi6 states that "^the construction of multitone signals having minimum peak-to- 
peak amplitude Is an old problem that still has no analytical solution". 

[0009] It should be noted that Golay complementary sequence pairs and codes defined using them are different from 
and should not be confused with so-called Golay codes; Golay complementary sequence pairs and Golay codes were 
first defined, independently of one another, by the same researcher and for that reason alone both are referred to by his 

30 name. 

[0010] Methods and apparatus tor encoding and decoding are described in international patent application no. 
PCT/GB 97/02 464 (WO 98/1 1 698), in which a substantial number of different code words are identifiable in practice, 
despite the potentially huge search space of possible code words, and which provide a predictable low limit on PMEPR, 
have specified error detection/correction properties and are feasible for implementation in high-speed circuitry. In par- 

35 tlcular codes are described for systems in wliich each transmitted symibol has 2^ different possible values where ; is a 
positive integer, such as binary 0' = 1). quaternary (/ = 2), octary (/=3). and higher orders. Techniques are also 
desaibed for decoding these codes, involving use of the fast Hadamard transform (FHT). These techniques are effec- 
tive, but they nonetheless entail substantial amounts of processing, and many applications of the f&sX Hadamard trans- 
form, especially in the case of octary and higher-order codes and longer code words (i.e. with more symbols). 

40 [001 1 ] It is an object of this invention to provide methods and apparatus for decoding such codes and which provide 
a reduction in the amount of processing required. 

Disclosure of Invention 

45 [001 2] According to one aspect of this invention there is provided a method of identifying a code word within a prede- 
termined set of code words, the code words containing symbols each having pi different possible values where p is a 
prime and j is an integer greater than 1 . wherein 

a plurality of constituents, constituent^ to constituent^, of said code word, each constituent identifying an element 
50 of a predetermined subset of a predetermined linear code, are derived in respect of an input word; and 
sakJ k constituents are combined to identify said code word according td a relationship 

(p^'^^*constituent;, + p''('''Vconstituent;t.i + ... + p^^^*constituent2 + p''^^^*constituenti) mod p^ 

55 where each b(x) is an integer, j > b[l<) > b{k-^) > ... > b{2) > 6(1) ^ 0. * indicates symbol-wise multiplication, + indi- 
cates symbol-wise addition and mod indicates symbol-vyise modular reduction. 

[0013] According to another aspect of this invention there is provided a method of identifying a code word within a 
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predetermined set of code words, the code words containing symbols each having 2^ different possible values where j 
is an integer greater than 1 . comprising the steps of 

receiving an input word; 

deriving the result of reducing the input word modulo 2^ ; 

deriving a first constituent, constituent^ , by applying a Hadamard transform operation, or equivalent thereof, in 

respect of the result of said reduction modulo 2^ ; 

setting a modifier equal to constituent^ ; 

for each successive value of a counter x from 2 to j: 



' deriving the resuK of symbol-wise subtraction of the current value of the modifier from the input word modulo 
2^: 

deriving another constituent, constituentj^ , by applying a Hadamard transform operation, or equivalent thereof, 
in respect of the result of said subtraction modulo 2^; 
15 - setting the modifier equal to the symbol-wise sum modulo 2^ of the modifier and the product of 2^'^ times 
constituent^;. 

and 

identifying the code word as having the final value of said modifier. 
20 [001 4] According to a further aspect of this invention there Is provided apparatus for identifying a code word within a 
predetermined set of code words, the code words containing symbols each having different possible values where p 
is a prime and j is an integer greater than 1 , cotrprising: 

a decoder for deriving in respect of an irput word a plurality of constituents, constituent^ to constituent/^, of said 
25 code word, each constituent identifying an element of a predetermined siA>set of a predetermined linear code; and 
a combiner for combining said k constituents to identify said code word according to a relationship 

{p^^'^^*constituent;f + p''('^"Vconstituent^.^ + ... + p^^^^constituentg + p''^^'*oonstituenti) mod p^p^ 

30 where each f?(x) is an integer.; > b[k) > b{k'^) > ... > b{2) > b(^) ^ 0. * indicates symtx)l-wise multiplication, + indi- 
cates symbol-wise addition and mod indicates symbol-wise modular reduction. 

[001 5] Decoders implementing this invention enable decoding to be accomplished with relatively few sqsplications of 
the fast Hadamard transform. For example, a quaternary code word can be decoded with just two successive applica- 
35 tions of the FHT, and an octary code word with just three, irrespective of the length of the code word. 

• * 

Brief Description of Drawings 

[0016] Methods and apparatus for decoding data in accordance with this Invention in communications systems for 
40 transmitting data using coded orthogonal frequency division multiplexing (COFDM) will now be described, by way of 
example, with reference to the accompanying drawings in whlch> 

Figure 1 illustrates multipath propagation between a base station and a mobile communications device; 

Figure 2 is a schematic diagram showing the effect of multipath propagation on propagation of a single sym- 

45 bol; 

Figure 3 is a schematic diagram showing intersymboi interference arising from multipath propagation; 

Figures 4a to 4h are schematic diagrams showing reduction in intersymboi interference by the use of multi-carrier 

frequency division multiplexing; 
Figure 5 shows an encoder and modulator for use in a COFDM system; 

so Figure 6 Is a block schematic diagrarn of a COFDM system incorporating the encoder and modulator of Fig- 

ure 5l 

Figure 7 Illustrates the general form of a COFDM signal waveform without any limitation of PMEPR; 

Figure 8 shows a decoder incorporating the present invention for decoding 16-symbol octary code words; 

and 

55 Figure 9 shows a decoder based on tiie decoder of Figure 8, for decoding 1 6-symbol octary code words and 

identifying an associated coset representative. 
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Best Mode for CarrvinQ Out the Invention. & Industrial Applicability 

[0017] Figures 1 to 3 illustrate in a schematic manner tiie problem of muttipath propagation in the context of mobile 
communications in an urban environment. A base station 1 0 communicates by wireless, e.g. at microwave wavelengths, 

5 with mobile communications devices such as a mobile telephone 12. There is a direct line-of-sight signal path 14 
between the base station 10 and the telephone 12, and, owing to the presence of multiple reflectors (for miaowave sig- 
nals) in tiie form of buildings and other structures, there are also multiple reflected signal paths 16. These paths are 
typically of differing lengths and thus involve differing propagation delays upon an'ival at their destination. The overall 
signal as received at. for example, the mobile telephone 12 is tiie sum of the various signals arriving at tiie telephone 

10 via tiie different patiis. 

[0018] The transmitted signal may be considered as comprising a succession of symbols modulating a carrier wave 
(e.g. by phase-shift keying). As shown in Figure 2, a consequence of signals traversing multiple paths with different 
propagation delays is that the composite signal 18 as received for a symbol is not an accurate replica of the signal 20 
as transmitted for that symbol. In particular tiie received composite signal is typically longer in duration than tiie trans- 

15 mitted signal, and has a trailing portion with a significant amplitude after the point when an accurate replica of tiie trans- 
mitted signal would have terminated. If this trailing portion extends for a significant fraction of a symbol period, as shown 
in Figure 3 for the pulse nominally ternninating at time T tiie ti-ailing portion will overlap and interfere with reception of 
the following symbol. As a result that following symbol may be incorrectly received, i.e. the value detected for it in the 
receiver may not correctiy represent tiie value as transmitted, even In the absence of other noise sources. Signal cor- 

20 ruption of tills kind is known as intersynribol interference (ISl). 

[001 9] The effects of ISl can be alleviated by extending the duration of the transmitted symbols, as shown in Figure 
4a, so that each symbol as received extends beyond the trailing portion arising from multipath propagation delays. To . 
counteract tiie reduction in symbol rate which would ensue, it has been proposed to transmit multiple symbols simulta- 
neously, each modulating a respective one of multiple carrier waves. Thus if eight carriers are used, as shown in Fig- 

25 ures 4a to 4/7, the symbol period on each carrier can be eight times tiiat which would be necessary with a single carrier 
(as shown in Figure 3) while maintaining the same overall symbol rate. If tiie frequency spacing between tiie multiple 
carriers is an integer multiple of the inverse of the symbol period; It is possible to recover the individual symbol stream 
modulating each of tiie can-iers witiiout interference from the other carriers. This technique is known as orthogonal fre- 
quency division multiplexing (OFDM). 

30 [0020] In addition to the use of multiple carriers, it has also been proposed to provide additional error detection/cor- 
rection capabilities by the use of block encoding. Thus, as shown in Figure 5, a serial stream of data symbols (binary 
symbols in this example) is notionally divided into successive groups or words each containing five symbols. A group of 
five symbols is input to a five-bit to eight-bit (5B/8B) encoder 30 which outputs a corresponding eight-bit code word. This 
code word is applied to a bark of eight phase shift controllers 32 which control a bark of eight osdilators 34 producing 

35 earner signals at eight successive frequencies at intervals Af apart. Each individual bit in the code word causes a 
respective phase shift oonta'oller to impose either a positive or negative phase shift A<p to modulate the carrier signal pro- 
duced by the respective oscillator, depending on tiie current value of that bit. 

[0021 ] The modulated outputs of tiie oscillators 34 are combined in an adder 36 to produce tiie transmitted signal. In 
practice the functions of the oscillators 34 and the adder 36 are conveniently combined in a digital signal processor 

40 inplementing an inverse fast Fourier ti^ansform (IFFT), as indicated by the dashed box 38. Thus, as shown schemati- 
cally in Figure 6, data to be transmitted aire received in block 40 and converted to parallel form in a serial -to-parallel con- 
verter 42. The parallel data are encoded in encoder 30. as described above, and the encoded data confrol generation 
of the multi-carrier signal in the modulator/1 FFT 32/38. This multi-carrier signal is converted to analogue form in a D-A 
converter 44, amplified and f ansmitted by a transmitter 46. 

45 [0022] Reception is essentially a reverse of tiiis process: the signal is received and amplified by a receiver 48. and 
then converted to digital form by an A-D converter 50. The digital signal is supplied to a demodulator 52 which com- 
prises a digital signal processor implementing the fast Fourier transform (FFT) to separate and denuxiutate the compo- 
nent signals in the composite multi-earner signal. The denruxiulated data are passed to a decoder 54 which reverses 
the encoding applied by tiie encoder 30, and typically also performs en'or detection or connection functions. Thereafter 

so a parallel-to-serial converter 56 changes tiie data into serial form to be output at block 58. 

[0023] The choice of code which is implemented by the encoder 30 has significant implications for tiie efficient oper- 
ation of tiie transmitter 46. A multi-carrier signal comprises tiie sum of several equal-amplitude sinusoidal waves having 
equally spaced frequencies. As a result, the maximum absolute value of the composite signal is sf ongly sensitive to 
the relative phases of tiie component sinusoids. These relative phases are changed at the end of each symbol period 

55 by the phase shift controllers 32. according to the value of the code word to be transmitted during the next symbol 
period. Depending on the particular set of relative phases, and therefore on the value of the current code word, ttie 
amplitude of tiie fransmitted signal will have occasional pronounced peaks together with intervals of much lower though 
still varying value, as illustrated in Figure 7. Thus the peak-to-mean envelope power ratio of the signal is relatively high. 
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[0024] For such a signal to be transmitted without distortion, rt is necessary for the transmitter 46 to use a linear ampli- 
fier and to operate it in a relatively inefficient manner, so that most of its dynamic range is unused except for the infre- 
quent occasions when the amplitude of the transmitted signal reaches its peak value. 

[0025] • It has been proposed to reduce this problem by appropriate choice of the set of code words which are actually 

5 used to encode data for transmission. The use of block coding to avoid the use of code words which produce extreme 
values of PMEPR has been suggested in 'Minimisation of the peak to mean envelope power ratio of multicanier trans- 
mission schemes by block coding' by T.A. Wilkinson and A.E. Jones. 1995 IEEE 45th Vehicular Technology 
Conference, July 1995, pp. 825-829; although this proposal Indicates that reductions in PMEPR are possible, it also 
identifies the difficulty of selecting a set of code words which both enables efficient inplemerrtation of coding and 

w decoding and provides desired capabilities for detecting/correcting errors in the demodulated signal. Possible ways of 
addressing this problem are discussed, but no actual solution is presented. A modification to this approach is disclosed 
in 'Combined coding for en-or control and increased robustness to system noniinearities In OFDM' by the same authors 
in 1996 IEEE 46th Vehicular Technology Conference, April 1996, pp. 904-908. In this modification a linear block code 
is chosen to provide desired properties of enror detection/correction; it is proposed then to exploit the linearity and 

15 redundancy in the code to transform the code values systematically to new values which reduce the PMEPR but have 
equivalent error detection/correction properties. Examples are given for 4BA7B. 4B/8B and 1 1 B/15B codes, in which the 
required transformation is identified by a procedure involving listing all possible code word values in order of the 
PMEPR for a multi-carrier signal modulated in accordance with each code word. Although such an exhaustive-search 
approach is feasible for relatively short code words (e.g. of the order of 15 bits), the number of possible code words to 

20 be checked increases very rapidly with code word length. Thus, even for a 32-bit code word. In excess of 4,000,000,000 
code words would potentially need to be listed, and this would be at best extremely time-consuming and in practice pos- 
sibly prohibitive. The obstacle presented by this requirement is potentially increased by the absence of any guarantee 
that an appropriate sut)set of code words actually exists within the set being examined. 

[0026] Popovi6 , in the paper dted earlier, proposes the use of a single binary or polyphase complementary sequence. 
25 However, the 1996 paper mentioned above suggests that such sequences are not readily amenable in relation to their 
error correction/detection capability. 

[0027] The patent application referenced above describes ways of identifying a set containing many different code 
words all of which satisfy desired criteria regarding error detection/correction properties, the maximum PMEPR of multi- 
can'ier signals modulated in accordance with those code words, and ease and efficiency of implementation of a practi- 
se cal encoder and decoder (for example including use of combinatorial logic). 

[0028] As described in that patent application, a set of code words is constructed by combining a 'base' code, derived 
from a 'generator matrix', with one or more 'coset repr^entatives'. For example, in the case of 16-bit (binary) code 
words the generator matrix is: 

^ (0000 0000 1111 1111) (Xi) 

(0000 11110000 1111) ' (X2) 

(0011 0011 0011 0011) (X3) 

« (0101010101010101) (X4) 

(1111 nil 1111 1111) (X') 

The 1 6-bit base code comprises all thirty-two linear combinations of rows of this generator matrix, that is combinations 
45 of the form 

(aiX, +a2X2 + a3X3-t-a4X4+'a'X^mod2 

where addition and modular reduction are performed bitwise and row coefficients ai, a2. as, 34 and a' each take values 
so of 0 or 1 . This base code is linear, i.e. the result of a bitwise exdusive-OR operation (modulo-2 addition) on any two of 
its code words is another code word in the code. 

[0029] A con-esponding generator array is created by combining the 16-bit generator matrix with additional rows (in 
this case six) derived by bitwise multiplication of all possible pairs of the rows X^, Xg, X3 and X4: 

55 
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(0000 0000 0000 nil) (X,*X2) 

(0000 0000 001 1 001 1) (Xi*X3) 

(0000 0000 0101 0101) (Xi*X4) 

(0000 001 1 0000 001 1) (X2*X3) 

(0000 0101 0000 0101) .(X2*X4) 

(0001 0001 0001 0001) (X3*X4) 



Twelve particular coset representatives (for example, 0001 0010 0001 1 101) constructed in a specific way from these 
additional rows are each combined by bitwise exclusive*OR operatior^ with all code words in the 16-bit base code to 
produce the final desired set of code words, constituting a set of cosets of the base code. 

[0030] Those code words constitute a set of cosets of a linear sub-code (in this case the whole code) of a base code 
having a generator matrix comprising the rows X^. X2, X3, X4 and X*. Using a base code of this Kind provides desired 
error detection/conrection properties. Selection of a linear sub-code ensures that encoding and decoding within each 
coset can be implemented compactly, for example using combinatorial logic. A coset can in general be considered to 
be a systematic shift or translation of the base oode to a different part of the overall space of possible code words. Use 
of cosets of the linear sub-code enables avoidance of certain undesired code words (such as all zeroes and all ones) 
which would result in particularly high values of PMEPR, but can also enable the error control properties of the base 
code to be retained to some extent. Choosing the code words from a set of such cosets allows good control of PMEPR 
because each coset of the base code tends to contain code words having broadly similar associated PMEPRs. In par- 
ticular cases where the coset representatives. have a form described in the patent application referenced above, it can 
be shown that the PMEPR does not exceed 3 dB; the minimum Hamming distance (i.e. the minimum number of sym- 
bols whose value must be altered, for example by interference, to transform a valid code word containing 2"^ binary 
symbols into some other valid code word) is 2"^'^. 

[0031] In more general terms the generator matrix for a particular code word length 2^ is: 

((2.0+l)x2'^/2'" ((2.1 + l)x2'^/2'" ((2.z-!-l)x2'^/2'" ((2.(2'"-l)+l)x2'^/2'" 
((2.0+l)x2'^')/2'^* ((2.1 + l)x2«'')/2« ' ((2.z+l)x2'^V2"*•^. ((2.(2"-l)+l)x2'"-V2'"-' 



((2.0+l)x20/2> ((2.H-l)x2W - ((2.z+l)x2W 



((2.(2'"-l) + l)x2W 



((2.0+l)x2V2* 
((2.0-l-l)x2°)/2° 



((2.1 + l)x2V2* 
((2.1-M)x2V2° 



((2.2+l)x2V2* 
((2.z4-l)x2°)/2" 



((2.(2'"-l)+l>x2V2' 
((2.(2«-l)+l)x2V2® 



where z has integer values from 0 to 2'"-1 , y has Integer values from m to 0, and x indicates bitwise multiplication (bit- 
wise logical AND) of two factors (such as (2,z+^) and 2^ both expressed as (m + 1)-bit binary numbers, and the result 
of division is expressed as a single digit. The first line of this matrix (y = m) gives the row X^ , the second line (y s m-1) 
gives the row X2, the penultimate line (y = 1) gives the row X^ and the final line (y = 0) gives the row X*. 
[0032] TTie 2'"-bit base code comprises all linear combinations of rows of the generator matrix, that is combinations 
of the form 

(a^X^ +a2X2 + 83X3 + .... + a^X^ +a'X') mod 2 

where row coefficients a^ , a2, sl^, • . a^^ and a' each take values of 0 or 1 . The corresponding generator array is created 
by combining the 2^-bit generator matrix with additional rows derived by bitwise multiplication of all possible pairs of 
the rows X^. X2, X3 .... X^. Particular coset representatives derived from these additional rows in a specific way are 
each combined by bitwise exdusive-OR operations with all code words in the 2'"-bit base code to produce the required 
code words, constituting a set of cosets of the base code. ^ 
[0033] "This technique is also applicable to coding with any higher number 2^ of symbol values, such as quaternary (j 
= 2; 2^ = 4), octary (j = 3; 2> = 8) and higher orders. In such cases the same generator matrix is used as described 
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above, e.g. for a code word containing 8 symbols {2"^ = 8 with m = 3): 



(0000 1111) (Xi) 

(00110011) (X2) 

(0101 0101) (X3) 

(1111 1111) (X') 



[0034] An 8-symbol quaternary base code (four possible symbol values) comprises all 256 linear combinations of the 
rows of this generator matrix, that is combinations of the form 

15 (a^X^ +33X2 +83X3 +a'X') mod4 

where addition is now performed symbol-wise modulo 2^ (i.e. modulo 4 for quaternary), modular reduction is performed 
symbol-wise and row coefficients a-i . a2. aa and a' each take all integer values from 0 to 2^-1 (i.e. 0 to 3 for quaternary). 
This base code is linear, i.e. the result of symbol-wise addition modulo 4 on any two of its code words is another code 
20 word In tiie code. The corresponding generator array is created by combining the generator matrix with additional rows. 
These additional rows are derived by symbol-wise mult^lication of all possible pairs of the rows X^, X2 and X3, that is 
X^ *X2. Xi *X3 and X2*X3, and tiien further multiplication by 2^*^ . e.g. by 2 for quaternary coding or Isy 4 for octary coding. 
Thus for 8-symbol quaternary coding the additional rows are: 



(0000 0022) 2(X,*X2) 
(0000 0202) 2(Xi*X3) 
(0002 0002) 2(X2*X3) 



[0035] Likewise a 2'"-symbol quaternary base code comprises all linear combinations of rows of the generator matrix, 
that is combinations of the form 

35 (a^X^ +a2X2 + a3X3 + .... + a^X^ + a'X7mod4 

where row coefficients a^, a2. blz ■■■ a^ and a* each take values of 0, 1 , 2 or 3. The corresponding generator array is 
created by combining ttie 2'"-symbol generator matrix with additional rows derived by symbol-wise multplication of all 

possible pairs of tiie rows X^, X2. X3 X^. Particular coset representatives derived from these additional rows in a 

40 specific way are each combined by symbol-wise addition modulo 4 with all code words in the 2'"-symbol base code to 
produce the required code words. 

[0036] Similar codes for otiier code word lengths and/or for coding with different numbers of symbols (e.g. octary) can 
be defined in like manner. Thus, where the number of symbol values is 2^, the base code comprises all linear combina- 
tions of rows of the generator matrix, tiiat is combinations of the form 



(a^X^ +32X3+33X3 + .... +a^X^ +aX)mod2^ 



where addition is performed symbol-wise modulo 2^ modular reduction is performed symbol-wise and row coefficients 
31 . 32, 33, .... a;„ and a' each take values of 0,1 2^-1 . This b3se code is linesr. i.e. the result of 3 symbol-wise addi- 

so tion modulo 2^ on sny two of its code words is snotiier code word in the code. The corresponding generstor array is cre- 
ated by combining the 2'"-symbol generator matrix with additional rows, which are derived by symbol-wise 
multiplication of ail possible pairs of the rows X^. X2, X3. .... X;„ and then further multiplication by 2>'^ Particulsr coset 
representatives derived from these addrtionsi rows in a specific way are each combined by symbol-wise addition mod- 
ulo 2^ with all code words in the 2^-symbol base code to produce the required code words. 

55 [0037] For a 1 6-symbol octary code (i.e. rn = 4 and j = 3), as discussed in examples given below, the generator an^ay 
is: 



Page 8 of 26 



8 



EP 924 910 



Page 9 of 26 



EP 0 924 910 A1 



(0000 0000 1111 1111) (X,) 

(0000 111100001111) (Xj) 

(0011001100110011) (Xj) 

(Old 0101 0101 0101) (X4) 

(1111 nil nil nil) (x') 

(0000 0000 0000 4444) 4(X,*X2) 

(0000 0000 0044 0044) 4(X,*X3) 

(0000 0000 0404 0404) 4(X,*X4) 

(0000 0044 0000 0044) 4(X2*X3) 



(0000 0404 0000 0404) . 4(X2*X4) 

(0004 0004 0004 0004) HX3*X^ 

and the twelve coset representatives constructed from these additional rows are: 

25 



Table 





(0004 0040 0004 4404) 


4(X^ *X2)©4(X2*X3)e4(X3*X4) 


30 


(0004 0400 0004 4044) 


4(Xi *X2)e4(X2*X4)©4(X3*X4) 


(0000 0440 0044 0404) 


4(Xi *X3)©4(X2*X3)©4(X2*X4) 




(0004 0400 0040 0444) 


4(Xi *X3)©4(X3*X4)e4(X2*X4) 




(0000 0440 0404 0044) 


4(Xi *X4)®4(X2*X4)e4(X2*X3) 


35 


(0004 0040 0400 0444) 


4{Xi *X4)®4(X3*X4)e4(X2*X3) 




(0004 0004 0040 4404) 


4(Xi *X2)e4(Xi *X3)©4(X3*X4) 




(0004 0004 0400 4044) 


4(Xi *X2)®4(Xi *X4)©4(X3*X4) 


40 


(0000 0044 0440 0404) 


4(X2*X3)®4(Xi*X3)©4(Xi*X4) 




(0000 0404 0440 0044) 


4(X2*X4)®4(Xi*X4)©4(Xi *X3) 




(0000 0404 0044 4004) 


4(Xi *X3)®4(Xi *X2)e4{X2*X4) 




(0000 0044 0404 4004) 


4(X2*X3)©4(Xi *X2)®4(Xi *X4) 



45 



Selected ones of these coset representatives (in the examples below, the first eight) are combined by symbol-wise addi- 
tion modulo 8 with all code words in the 1 6-symbol base code to produce the final set of octary code words. 
[0038] In discussing the error detecting/correcting properties of quaternary, octary and higher-order codes, the mini- 

50 mum Hamming distance alone does not filly describe the extent of possit)le errors affecting a code word. This is 
because an error in a single synlboi can involve a transition through up to 2^'** adjacent symk>ol values (where the values 
2^-1 and 0 are considered to be adjacent). The practical effect may be as serious as. or more serious than, two or more 
errors in different symbols each involving only a single transition to an adjacent symbol value. For example, a change 
of a single octary symbol from value 6 to vjalue 1 (a transition through 3 adjacent values) may be more serious than a 

55 change of two octary symbols, one from 6 to 7 and another from 7 to 0 (both being transitions to an adjacent value). 
Therefore an additional useful measure is the minimum Lee distance, defined as the minimum number of symbol value 
transitions required to transform a valid code word into some other valid code word. In the binary case the Hamming 
and Lee distances are identical. 
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[0039] The methods and apparatus tor encoding and decoding COFDM signals as described in the above-referenced 
patent application provide a highly advantageous combination of error detection/correction properties, low PMEPR. 
convenience of implementation and useful rate of data transfer. Both encoding and decoding can be accomplished 
using circuitry which incorporates combinatorial logic, and the above-referenced patent application describes several 
5 such decoders based on Lee distance. 

[0040] These decoders are based on the "Supercede" decoding method described in 'Soft decoding techniques for 
codes and lattices, including the Golay code and the Leech lattice* by J.H. Conway and NJ.A. Sloane, IEEE Trans, 
Infom. Theory, 1986, vol.lT-32, no.1 , pp. 41 -50, but with a significant modification to enable its use with quaternary (and 
higher-order) codes. For each coset representative in turn that coset representative is subtracted symbol-wise (modulo 
10 2^ from the received code word, and the result is decoded as an element of the relevant base code, at the same time 
assigning a "score" to that coset representative. The decision as to which code word to decode to is made by selecting 
a coset whose representative attains the highest score (arbitrarily If more than one attains that score), and within that 
coset selecting the element of the base code previously identified for that coset. 

[0041] The decoding of an element of the base code uses the fast Hadamard transform (FHT), as described in The 
15 Theory of Error-Correcting Codes' by FJ. MacWilllams and N.J. A. Sloane, North-Holland. Amsterdam, 1986. Specifi- 
cally, the fast Hadamard transform Is applied to multiple input vectors whose symbols are not constrained to have the 
value + 1 or -1 (as is conventionally the case), in order collectively to decode the result of the above subtraction as if it 
were a received code word of the base code. For each input vector a transform coefficient of maximum magnitude is 
selected. According to the sign of that transform coefficient, one or other of an associated pair of elements of the first- 
20 order binary Reed-Muller code RM(1 .m) could then be output. Although this is analogous to the procedure for decoding 
a received code word of RM(1 .m), the output for each input vector would not represent a decoding of that vector as an 
element of RM(1./77). 

[0042] Although these decoders are effective, for higher-order codes (e.g. octary) and longer code words, the amount 
of processing involved, including the number of required applications of the FHT. could become prohibitive in practical 
25 situations. 

[0043] Figure 8 shows an alternative decoder 60. in this example for use with l6-symbot octary base code words, 
which significantly reduces the required number of applications of the FHT The individual blocks of this decoder 60 may 
In turn be incorporated in a decoder 70. shown in Rgure 9, which provides an advantageous variation of the modified 
Supercede method described above to decode a code word which has been combined with a coset representative prior 

30 to transmission. In particular, the decoder 70 interleaves the subtraction of coset representatives with the operation of 
the constituent blocks of the decoder 60, to give a much faster decoder than wouki be obtained by using the full proce- 
dure implemented in the decoder 60 in respect of each coset representative in turn. A score is assigned to each coset 
representative and used to select a coset representative whose score either attains the highest value, or exceeds a pre- 
determined threshold (in which case processing of further coset representatives is omitted). The decoders 60 and 70 

35 also apply the FHT to input vectors whose symbols are not constrained to have the value + 1 or -1 ; however, the FHT 
is used for a different purpose from that described in the above-referenced patent application. 
[0044] Referring to Figure 8, a 1 6-symbol octary received code word Tq to be decoded by the decoder 60 is provided 
to a first FHT circuit 62, which applies the FHT in respect of that word and derives, as described below, a 16-component 
vector fo (which is a constituent of the output code word) and a modified word r^ . 

40 [0045] The modified word r-i is supplied to a second FHT circuit 64, which applies the FHT in respect of the word r-| 
and derives a 16-component vector f ^ (another constituent of the output code word) and another modified word r2. This 
modified word r2 is supplied to a third FHT circuit 66. which applies the FHT in respect of that word and derives a third 
1 6-component vector f2 (a third constituent of the output code word). 

[0046] Each FHT circuit 62, 64 and 66 includes a so-called 'Green machine' for performing the FHT (see the book by 
45 MacWilliams and Sloane referenced above, pp. 423-5). or a modification of the FHT which identifies a transform com- 
ponent of largest absolute value without deriving all transform components (see. for example, 'Fast decoding algorithms 
for first order Reed-Muller and related codes* by A.E. Ashikhmin and S.N. Utsyh, Designs, Codes and Cryptography, 
vol.7. 1996. pp. 187-214). Alternatively the FHT circuits 62, 64 and 66 may identify a transform component of largest 
absolute value by calculating a Hadamard transform (or a modification analogous to that mentioned above) without 
so using the actual FHT algorithm itself. 

[0047] The three vectors fo. f i and f2 are supplied by the FHT circuits 62, 64 and 66 to a combiner circuit 68 which 
combines the constituent vectors to derive a decoded 1 6-symbol output word F according to the relationship 

F = (4f2 + 2fi + fo) hiodS 

55 

[0048] The operation of the circuits 62 to 68 in the decoder 60 will be desaibed in more detail by reference to a spe- 
cific example of a 1 6-symbol octary base code word 
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(a ^ X ^ + a 2X 2 + 83X3 -t- 84X4 -f a'X*) mod 8 

with row coefficients a1 = 5, 83 = 7, 83 = 3, 84 = 6 and a' = 6, i.e.with a value of 6417 5306 3164 2053. This code word 
is assumed to have been corrupted in transit and has been received as 241 1 6316 3152 2453, which constitutes the 
5 value for the octary word Tq. 

[0049] In the first FHT circuit 62 a 1 6-connponent vector Dq is derived by applying the expression 

|1-(ro)z mod 21-1/2 (1) 

10 for values of z from 0 to 15 to define each respective component of the vector Dq. In this expression (1) the term (ro); 
indicates conrponent z of the word Tq, It will be noted that expression (1) yields values of +1/2 or -1/2; for convenience 
of implementation, each component may In practice be multiplied by two without affecting the derived values fo and r-i. 
For the example (enrored) received code word given above, these components of Dq (after multiplication by two) are 

15 (1, 1. -1. -1. 1. -1, -1. 1, -1, -1. -1. 1. 1. 1. -1. -1) 

[0050] The first FHT circuit 62 then applies the FHT to the vector Dq to derive sixteen transform components labelled 
from 0 to 15. as follows: 

20 (-2. -2, 6, 6. -2. -2, -2. -2. 2. 2. 2. 2. 2. 2. 10,-6) 

A component of largest absolute value (in this case, 10) is identified, the label (14) for this component is expressed as 
a four-bit binary value C1C2C3C4 (1110), and a coefficient c' is set to a value of 0 or 1 according to whether the component 
Is positive or negative respectively (in this case, 0). Then the vector fo is derived according to the relationship 



25 



55 



f 0 = (c iX 1 + C2X2 + C3X3 + C4X4 + c'X') mod 8 



which for the example received code word is (X ^ + + X3) mod 8 = 001 1 1 122 1 122 2233 . Here and in all subse- 
quent instances of addition and/or modular reduction of vectors, these operations are performed symbol-wise. 
30 (0051 ] The modified word r-i supplied to the second FHT circuit 64 is derived according to the relationship 

where the subtraction of vectors is performed symbol-wise (here and subsequently). Thus in this case r-i has a value of 
35 2400 5274 2030 0220. 

[0052] In the second FHT circuit 64 a 1 6-component vector is derived by applying the expression 

|2-(ri)^mod4|-1 (2) 

40 for values of z from 0 to 15 to define each respective component of the vector In this expression (2) the term (r^)^ 
indicates component z of the word r^. For the current exanple, these components of are 

(-1.1.1.1,0. -1,0,1.-1.1.0,1.1.-1,-1.1) 

45 The second FHT circuit 64 then applies the FHT to the vector to derive sixteen transform components labelled from 
Oto 15, as follows: 

(3, -5, -5, 3. 3. -5. -1, -9, 1, 1. -3. -3, 1. 1. 1. 1) 

50 A component of largest absolute value (-9) is identified, the label (7) for this component is expressed as a four-bit binary 
value C1C2C3C4 (01 1 1), and a coefficient c* is set to a value of 0 or 1 according to whether the component is positive or 
negative respectively (in this case, 1). Then the vector f^ is derived according to the relationship 



f^ = (c^Xi t C2X2 + C3X3-i-C4X4 + cX)nrK)d4 

i.e. (X2 -f X3 + X4 + X*) mod 4 = 1223 2330 1223 2330 . The modified word rg supplied to the third FHT circuit 66 is 
derived according to the relationship 
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r2B(r^ -2f^)mod8 

Thus in the present example r2 = 0042 1414 0472 4440. 

[0053] In the third FHT circuit 66 a 16-component vector D2 is derived by applying the expression 

|4-(r2)^mod8|-2 (3) 

for values of z from 0 to 15 to define each respective component of the vector D2. In this expression (3) the term (rg); 
indicates component z of the word rg. For the present example these components of D2 are 

(2, 2, -2, 0. 1, -2, 1. -2. 2. -2, 1. 0. -2, -2. -2. 2) 

[0054] The third FHT circuit 66 then applies the FHT to the vector D2 to derive sixteen transform components labelled 
from 0 to 15, as follows: 

(-3, 5. 1,9, 9. 1.9.1.3. 3. 11. -5, -1.-17. 3.3) 



A component of largest absolute value (-17) is identified, the label (13) for this component is expressed as a four-bit 
binary value 0^020304 (1101). and a coefficient C is set to a value of 0 or 1 according to whether the component is pos- 
20 itive or negative respectively (in this case. 1). Then the vector f2 is derived according to the relationship 

fg = (CiXi +C2X2 + C3X3 +C4X4 +c*XO mod 2 

i.e. (Xi +X2 + X4 + X') mod 2= 1010 0101 0101 1010. 
25 [0055] The three vectors fo, f 1 and are supplied to the combiner circuit 68 for further processing as described above 

to derive an output code word F For the present example F o (4(1010 0101 0101 1010) + 2(1223 2330 1223 2330) + 

(001 1 1 122 1 1 22 2233)) mod 8 = 6417 5306 3164 2053, I.e. the base code word as transmitted. 

[0056] If tt is desired to recover the row coefficients a^ . a2. 33. 84 and a' then F is evaluated In terms of the generator 

matrix rows X^ . X2. X3, X4 and X* rather than as a 16-component vector. Thus, in the above example, F would be eval- 
30 uated as (4(Xi + X2 + X4 + X') + 2(X2 + X3 + X4 + X*) + (Xi + X2 + X3)) mod 8 = {5X ^ + 7X g + 3X 3 + 6X 4 + SX") mod 8 . 

In this case the procedure stili identifies the code word F even though F need not be explicitly determined as part of the 

procedure. 

[0057] The procedure described with reference to Figure 8 can be considered as identifying elements, fo, f-i and fg. 
each belonging to a respective predetermined subset of a predetermined linear code. In this procedure every element 
35 of such a subset is a candidate result for fo, fi or fg as appropriate. In the case of fo, for example, the predetermined 
linear code is the 16-symtx)i octary base code defined above, and the predetermined subset is the set of all vectors of 
the form 



(c^X., + C2X2 + C3X3 +C4X4 + c'X') mod 8 

in which each of the coefficients Ci , eg, C3. C4 and c* can take only the values 0 or 1 (rather than the values 0 to 7 inclu- 
sive). For f^. the predetermined linear code Is the 16-symbol quaternary base code defined above, and the predeter- 
mined subset is the set of all vectors of the form 

45 (c^X^ +C2X2 + C3X3 + C4X4+c'X')mod4 

in which each of the coefficients c^ , C2. C3. C4 and c* can take only the values 0 or 1 (rather than the values 0 to 3 inclu- 
sive). 

[0058] The decoding method embodied in the decoder of Figure 8 can readily be extended for use with other lengths 
so of code word and/or higher-order codes. In the general case of code. words of length 2^" symbols, each symbol having 
2^ possible values, a total of j FHT circuits are required, which for the purposes of. the following desaiption are identified 
by an index k having values from 0 (con-esponding to the first FHT circuit 62 in Figure 8) to y-1 (corresponding to the 
third FHT circuit 66 in the octary case). 

[0059] In the FHT circuit with index k a vector D/^ is derived from the word tf^ received by that circuit, by applying the 
55 expression 

|2''-(r/,)^mod2''*V2''*^ (4) 
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for values of z from 0 to 2'"-1 to d^ine each respective component of the vector D^. In this expression (4) the term (rj^)^ 
indicates component z of the word r^^. 

[0060] Each FHT circuit then applies the FHT to the vector D to derive 2^" transform components, labelled from 0 to 
2!^ A . A component of largest absolute value is identified, the label for this component is expressed as an m-bit binary 
5 value C1C2... C;i„ and a coefficient C Is set to a value of 0 orl according to whether the con^X)nent is positive or negative 
respectively. Then the vector f is derived according to the relationship 

f^ = (CiXi +C2X2 + ... + c^X^+c'X')mod2^'^ 

10 For all but the final FHT circuit [k = /-I) the modified word r^^^ supplied to the next FHT circuit is derived according to 
the relationship 

r/r+i=(r;^-2'^f;c)mod2^ 

IS [0061 ] The y vectors fo. f 1 . • f/.i are combined in a combiner circuit (corresponding to the circuit 68 for the octary case) 
according to the relationship 

F = (2^'^fy.i + 2^'^fy.2 + ... + f 0) mod 2^ 

20 to produce the required output base code word F. 

[0062] it should be noted that in view of the repetitive pattern of the derivation in each FHT circuit, distinguished solely 
by the value of the index k from circuit to circuit, the decoding method may also readily be implemented in an Iterative 
manner, for example by a processor executing appropriate software program, instructions. 
[0063] Furthermore, the relationship for the output base code word F can alternatively be expressed as 

25 

F = (2{...2(2(2(2fy.i + fy.2) + f y.3) + fy^)... + f 1) + f 0 mod 2^ 
Thus F can clearly be evaluated in a recursive manner. For example, in the octary case discussed above 

■ 

30 F = (2(2f 2 + f 1) + f 0) mod 8 

The vector fo is derived as described above, using the first FHT circuit 62. Then the output code word F is considered 
to be the symbol-wise sum modulo 8 of fo and twice tiie output code word of a decoder for tiie quaternary case, in place 
of the FHT circuits 64 and 66. Likewise, for the case j « 4, the output code word is considered to be tiie symbol-wise 

35 sum modulo 16 of fo and twice tiie output code word of a decoder for the octary case, which in turn is treated either 
recursively as just described, or directly as described earlier. Other modifications of a similar nature can be envisaged. . 
[0064] The decoder as described above is a 1iard decision' decoder, i.e. it takes integer values between 0 and 2^-1 
inclusive as input. The same decoder can also be used for 'soft decision' decoding, wherein tiie symbols of the received 
code word are represented as real numbers between 0 and 2^ according to the value of the received analogue signal. 

40 In deriving vectors 0;^ by evaluating expression (4) for real-valued elements {r/f)^, the term (r/^)^ mod 2^*'^ is taken to be 
the real number a having a value such that 0 ^ a < 2^^^ and such that (r^^)^ - a is an integer multiple of 2^^"^^ Each FHT 
circuit produces real-valued transform components, of which one with largest absolute value is selected. The rest of tiie 
decoding procedure is tiie same as described above. 

[0065] For decoding code words comprising a base code word combined with a coset representative, as described 
45 above, the decoder 70 shown in Figure 9 may be used. This decoder is intended for use witii 16-synibo! octary code 
words which comprise a code word selected from the octary base code and added symbol-wise modulo 8 to a selected 
one of tiie first eight coset representatives in Table 1 above. Encoders for implementing these steps are described In 
the above-referenced patent application. 

[0066] As noted above, that patent application also describes decoders which use a nrKXtifled version of the Super- 
so code decoding metiiod, in which each possible coset representative is 8ut>tracted symbol-wise (nrKXiulo 2^) from ttie 
received code word, and tiie result is decoded as an element of the base code, at tiie same time assigning a '^score" 
(tiie largest absolute value of the FHT components) to tiiat coset representative. Witii tiie decoding method illustrated 
in Figure 8 and coset representatives which consist solely of symbol values of 0 and 4, tiiere is no need to repeat tiie 
steps performed in the first and second FHT circuits 62 and 64 for each possible coset representative. These circuits 
55 derive the vectors Dq and D-i by modulo 2 and modulo 4 arithmetic respectively, and any coset representative compris- 
ing only symbol values of 0 and 4 reduces to a value of all 0 symbols under these arithmetics. Thus the subtraction of 
each coset representative need only be effected in respect of tiie third FHT circuit 66; the inputs to the first and second 
FHT circuits 62 and 64 need not involve tiie coset representatives. 
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[0067] Furthermore, It is not always necessary to repeat the processing in the third FHT circuit 66 for every relevant 
coset representative, with a view to comparing the scores and selecting a coset representative which attains the highest 
score. It is possible to determine both a maximum attainable score (maximum absolute transform component value), 
and a threshold which if exceeded for a particular coset representative causes the decoder 70 to onnit processing of fur- 
5 ther coset representatives. 

[0068] In the case of coset representatives derived from the generator array as described above, such as the coset 
representatives in Table 1 , this threshold is equal to 3/4 of the maximum absolute transform component value for the 
relevant FHT circuit. The maximum value and associated threshold for each drcurt are independent of the value of j and 
are as follows: 

10 



15 



Circuit 


Max 


Threshold 


FHT circuit 64 
FHT circuit 66 
FHT circuit with index k 


16 

32 


12 
24 

2 ^fn^k-^ 



20 [0069] The decoder 70 shown in Figure 9 takes advantage of these properties to restrict the amount of processing 
required. 

[0070] Referring to Figure 9 the 1 6-symbol odary code word received by tiie decoder 70 constitutes the input word 
ro for the first FHT circuit 62. As described above, tills circuit applies the FHT and derives tiie vector fo and a modified 
word Tv This modified word r-i is supplied to tiie second FHT circuit 64, which derives the vector f^ and tiie modified 
25 word r2. 

[0071] The modified word r2 is supplied to a moduIo-8 subtracter 72, which subtracts from that word, symbol-wise 
modulo eight, a value in a coset store 74. This coset store contains the values of the eight coset representatives from 
which selection is made during encoding of the code word. The store 74 also has a three-bit counter 76 which cycles 
through all eight possible three-bit values successively In response to an increment signal at an Input 78, so tiiat each 
30 of the eight coset representatives in turn is provided to the subtracter 72. 

[0072] TTie word output by the subtracter 72 is supplied to the third FHT circuit 66. which derives tiie third vector i2- 
The output of tiie Green machine in tiiis third FHT circuit 66 is a vector E, comprising sixteen transform components Eq 
to El 5. 

[0073] The vectors fg and E are supplied to a comparator 80 which compares the absolute values of the transform 
35 connponents in the vector E, and identifies the largest absolute value jE^ J; if more than one component has the same 
largest absolute value, one of them is chosen aifoitrarlly. The comparator 80 then compares the absolute value of the 
chosen component witii tiie appropriate tiireshold value described above, in this case a value of 24; if the tiireshold Is 
exceeded, no furtiier coset representatives need to be tested. In either case, the comparator 80 outputs the vector f2 
and the value jE^I to a latch 82. which also receives the corresponding tiiree-bit index value supplied by tiie counter 76 
40 in the store 74. 

[0074] For the first coset representative the latch 82 simply stores the vector f2, the value |E^| and tiie coset index 
which it receives. If the threshold has not been exceeded the counter 76 is incremented, and another vector f2 and value 
lE^I are derived and provided to the latch 82 for the next coset representative, togetiier with the respective coset index. 
If such a transform component absolute value [E^l is larger tiian the absolute value currentiy latched, tiie latch 82 

45 selects the newly received absolute value, together with its associated vector f2 and the coset index. This cycle is 
repeated for each successive coset representative, until eitiier tiie tiireshold Is exceeded or the counter 76 has cycled 
through all eight possible values and- every coset representative in tiie store 74 has been supplied to tiie subbBctor 72. 
[0075] Thereafter the two vectors fo and f ^ derived by the FHT circuits 62 and 64. and tiie vector f2 latched by the latch 
82. are supplied to tiie combiner circuit 68 to derive the output 16-symboi code word F from the base code according 

so to tiie relationship 

F = (4f2 + 2f^ +fo) mods 

The latched coset index identifies the coset representative which was combined witfi this base code word F in the orig- 
55 inai transmitted code word, provided the code word as received has not sustained excessive corruption. 

[0076] The operation of the decoder 70 will be illustrated by reference to the same example 16-symbol octary code 
word used above (6417 5306 3164 2053), which is assumed to be combined witii the tiiird coset representative in the 
store 74(0000 0440 0044 0404) to produce a transmitted code word 641 7 5746 31 20 2457. This transmitted code word 
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is assumed to have been corrupted in transit and has been received as 6437 5046 4320 3457 (=ro). 

[0077] In the first FHT circuit 62 the 16-component vector Dq is derived as described above, with the following com- 
ponents (after multiplication by two) 

The first FKU circuit 62 applies the FHT to this vector Dq to derive sixteen transform componerrts labelled from 0 to 15, 
as follows: 

10 (2.-2. 2.-2, 2. 6. 2. 6. 2.-2. 2.-2.-6.-2.1 0.-2) 

The component of largest absolute value (10) is identified, the label (14) for this component is expressed as a four-bit 
binary value c^c^c^ (1 110), and the coefficient c' is set to a value of 0. Then the vector fo is derived according to the 
relationship 



15 



20 



25 



30 



35 



SO 



io-(c-^X^ +C2X2-»-C3X3 + C4X4 4-c'X')mod8 

which in this example is (X^ + X2 + X3) mod 8 = 0011 1122 1122 2233. The modified word r-i supplied to the second 
FHT circuit 82 is derived according to the relationship 

Ti =(ro-^o)"^od8 

and has a value of 6426 4724 3206 1 224. 

[0078] In the second FHT circuit 64 the 1 6-component vector D-i is d^ived with components 

(-1.1.-1.-1.1.0.-1.1.0.-1.1.-1.0.-1.-1.1) 

The second FHT circuit 64 applies the FHT to the vector to derive sixteen transform conriponents labelled from 0 to 
15, as follows: 

{-3.:1. 1, 3.-3. 3. 1,-9. 1,-5. 5.-1,-3,-5. 1,-1) 

A component of largest absolute value (-9) is identified, the label (7) for this component is expressed as a four-bit binary 
value C1C2C3C4 (0111). and the coefficient c' is set to a value of 1 . The vector f^ is derived according to the relationship 



i\ =(0^X1 +C2X2+C3X3 + C4X4+c'X')rnod4 
i.e. (X 2 + X 3 + X 4 + )C) mod 4 = 1 223 2330 1 223 2330 . The modified word r2 is derived according to the relationship 
40 r2 = (ri -2fi)mod8 

and has a value of 4060 0 1 44 1 640 5444. 

[0079] The subtracter 72 subtracts the first coset representative 0004 0040 0004 4404 in the store 74 from the mod- 
ified word r2, to derive a word 4064 0104 1 644 1 040, which is supplied to the third FHT circuit 66. 
45 [0080] This third FHT circuit derives the vector D2 with components 

(-2, 2, 0,-2, 2. 1. 2.-2, 1. 0,-2.-2, 1, 2,-2. 2) 

and then applies the FHT to the vector D2 to derive the sixteen transform oorrponents Eq to E^s, as follows: 

( 1,-1.13,-5.-11,-1, 1.-5, 1, 7.-3.-13. 1.-13.-3,-1) 



A componerrt of largest absolute value |E^| (13) is identified, the label z (2) for this conrponent is expressed as a four- 
bit binary value C1C2C3C4 (0010), and the coefficient c'' is set' to a value of 0. The vector f2 is derived according to the 
55 relationship 

f2 = (c:|X^ +C2X2 +C3X3 + C4X4 +c'X') mod 2 
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i.e. (X3) mod 2 ^ 0011 001 1 0011 0011 . The latch 82 stores this vector ^2. the value 13 and the coset index 000. 

[0081] As the value jE^j has not exceeded the threshold value of 24. the next coset representative 0004 0400 0004 
4044 is supplied by the store 74 to the subtractor 72, which derives a word 4064 0544 1 644 1 400 and supplies it to the 
third FHT circuit 66. This circuit derives the vector D2 with components 

5 

(-2. 2. 0.-2. 2.-1.-2.-2. 1. 0.-2.-2, 1.-2, 2. 2) 

and applies the FHT to this vector to derive the sixteen transform components Eq to E^s, as follows: 

10 (-5. 5. 7, 1.-5.-7. 7.-11,'5,-3. 7,-7. 7.-3.-13.-7) 

A component of largest absolute value jE^j (-13) is identified, the label z (1 4) for this component is expressed as a four- 
bit binary value C1C2C3C4 (1 110). and the coefficient c' is set Id a value of 1 . The vector fg is derived according to the 
relationship 

15 

f 2 = (c + C2X2 + C3X3 + C4X4 + c'X^ mod 2 

i.e. (X 1 + X2 + X3 + X") mod 2 » 1 100 001 1 001 1 1 100 . As the value jE^^I (13) is not larger than the currently latched 
value (1 3), the contents of the latch 82 are not changed. 
20 [0082] Again the value |E^| has not exceeded the threshold value of 24, so the next coset representative 0000 0440 
0044 0404 is supplied by the store 74 to the subtractor 72, which derives a word 4060 0504 1604 5040 and supplies it 
to the third FHT circuit 66. This circuit derives the vector D2 with components 

(-2. 2. 0. 2. 2.-1. 2.-2. 1. 0, 2.-2.-1. 2.-2, 2) 

25 

and applies the FHT to this vector to derive the sixteen transform components Eq to £^5. as follows: 

(5,-1 , 1 ,-5. 1 ,-1 ,-3.-5. 1 . 3,-3.-1 . 1 ,-25.-3. 3) 

30 A component of largest absolute value lE^I (-25) is Identified, the label z (1 3) for this component is expressed as a four- 
bit binary value C1C2C3C4 (1 101). and the coefficient c' is set to a value of 1 . The vector f2 is derived according to the 
relationship 

fg = (c^X^ +C2X2 +C3X3 + C4X4 +c*X') mod2 

35 

i.e. (X^ + X2 + X4 + XT mod 2 = 1010 0101 0101 1010 . Because the value of |E^| (25) is larger than the currently 
latched value (13), the latch 82 stores this vector f2. the value 25 and the coset index 010. 

[0083] This time the threshold value of 24 has been exceeded, so it is not necessary to derive values of f2 for any 
further coset representatives. The three vectors fot f^ and f2 are supplied to the combiner circuit 68 as described above 
40 to derive the code word F = (4(1010 0101 0101 1010) + 2(1223 2330 1223 2330) + (001 1 1 122 1122 2233)) mod 8 = 
641 7 5306 3164 2053, i.e. the base code word incorporated in the original, unconupted code word as transmitted. This 
code word F and the coset index 01 0 comprise the output of the decoder 70. 

[0084] In the general case of code words of length 2^ symbols, each symbol having 2^ possible values, the subtractor 
72 subtracts each coset representative from the word ry.^ modulo 

45 [0085] As described in the above-referenced patent application, the coset representatives may be chosen to include 
symbols with values other than 0 and 2^'^ ; or example, for octary code words coset representatives which may include 
synfibol values of 2 and 6 as well as 0 and 4 may be chosen, to obtain a larger family of code words (at the expense of 
increased PMEPR). In such a case additional subtraction is needed before the second FHT circuit 64, with an additional 
comparator and latch after it analogous to the comparator 80 and latch 82. The complete set of chosen coset represent- 

50 atives is first reduced modulo 4, and all the distinct resulting residues are identified. It is these distinct residues which 
are subtracted, symbol-wise modulo 4. in turn in the additional subtraction operation to produce the input in respect of 
the second FHT circuit 64. The additional comparator and latch after this FHT circuit operate to select one of these res- 
idues, either because a score assigned to it attains the highest value, or because it exceeds a predetermined threshold, 
derived as described above. Processing in sut>sequent FHT circuits in the decoder continues in respect only of the 

55 coset representatives whose symbol-wise reduction modulo 4 equals the selected residue. 

[0086] The decoders described above may be implemented in special-purpose hardware; in this case, whenever two 
or more distinct residues of tiie coset representatives are to be processed (for exanple. in the third FHT circuit 66 in the 
decoder 70 of Figure 9). some or all of the associated fast Hadamard transforms can be calculated in parallel, for exam- 
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pie by providing a Green machine for each one. 

[0087] Modified versions of the decoders 60 and 70 desaibed above can correct certain en'or patterns which cannot 
always be corrected by the decoders 60 and 70 themselves (and vice versa). In these modified versions the (/c + 1)th 
FHT circuit calculates the vector in the following manner. A vector v;^ is derived from the word r^^ by applying the 
5 expression 

((r^)^mod2'^*')/2'^ 

for values of z from 0 to 2'"-1 to define each respective component of the vector v/^. The vector D;^ is then derived from 
70 the vector v/^ by applying the expression 1 -2(v/f); for values of z from 0 to 2'"-1 . 

[0088] Next each FHT circuit applies the FHT to the vector D/f to derive 2'" transform components, labelled from 0 to 
2^ A . A component of largest at^lute value is identified, the label for this component is expressed as an m-bit binary 
value C1C2.... Cm. and a coefficient c' is set to a value of 0 orl according to whether the component is positive or negative 
respectively Then the vector f^ is derived according to the relationship 

15 

f^ = (CiXi +C2X2 + ... + CmXm + c*X')mod2^"'^ 

For all but the final FHT circuit {k s >1) the modified word r^^^^ supplied to tiie next FHT circuit is derived according to 
the relationship 

so 

where the vector e/^ is derived as 
25 (v;^ + fjif) mod2 

[0089] The j vectors fo. fi . ... fy.i are combined in the combiner circuit 68 according to the relationship 

^ = (2^'^*y.i + 2^'^f y.2 + ... + f 0) mod 2^ 

30 

to produce the required output base code word F 

[0090] When using such a modified decoder in conjunction with coset representatives including symbol values other 
than 0 and 2^'^ (for example, symbol values of 2 and 6 as well as 0 and 4 for octary code words) the distinct residues 
modulo 2^*'^ of the remaining coset representatives are each subtracted from the vector r/^ prior to derivation of the vec- 

35 tor Ok in the (/c + 1)th FHT circuit. 

[0(^1] The discussions above of octary codes have been expressed in terms of octary symt}ols having possible val- 
ues of 0 to 7. For actual transmission, for example by 8-phase^hift keying (8-PSK) modulation, these values would cor- 
respond to eight different values of phase shift. One possibility is phase shifts of 0°. 45°. 90"*, 135*", 180". 225**. 270** 
and 315° (equivalent to complex vectors 1, (I+O/V2, /, (-1+/)/V2, -1, (-1-/)/V2, -/and {1-/)/V2 respectively where i=:V-1). 

40 However, other variants are possible in which the phase shifts are, for example, 22^A°, 67^/i°, 1 12V4**, 157!/^°, 2021/^°, 
247W, 292^y4° and 337^/^° (equivalent to complex vectors [V(2+>fe)+W(2-V2)]^. [V(2W2)+/V(2+V2)]/2. [-V(2- 
V2)+/S/(2+V2)]/2, [W(2+V2)+/V(2W2)]^, [W{2+V2)-/V(2W2)]/2, [W(2W2)-W(2+>/2)]/2. [V(2-V2)-/V(2+V2)]/2 and [V(2+^)- 
/V(2-V2)]/2, respectively). This principle applies for codes with more symbol values (/>3), and indeed in other applica- 
tions the symbol values may be represented by any 2^ distinguishable signal modulations; tiieir amplitudes need not be 

45 equal, and, in tiie case of phase-shift keying, the phase differences between adjacent pairs need not all be equal. 
[0092] The description above has referred to linear 2^-ary codes in which code word symbols each have 2^ possible 
values, such as quaternary (2^) and octary (2^), and the decoding operations are performed with respect to a 2-ary 
(binary) code. However, the invention is also applicable to p'-ary codes in which tiie code word symbols each have pf 
possible values, where p is a prime other than 2, for example a 3^-ary (nonary) system, and the decoding operations 

so are performed with respect to a p-ary code. 

Claims 

1 . A method of identifying a code word within a predetermined set of code words, the code words containing symbols 
55 each having pi different possible values where p is a prime and j is an integer greater than 1. wherein 

a plurality of constituents, constituent^ to constituent;^, of said code word, each constituent identifying an ele- 
ment of a predetermined subset of a predeterrnined linear code, are derived in respect of an input word; and 
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said k constituents are combined to identify said code word according to a relationship 

(p^^'^^*constituent;f + p^^'^'^^*constituent ^^.^ + ... + p^^^^conslituenta + p^^^^*constituent i) mod p^ 

where each b{x) is an integer, j > b(k) > /7(/c-1) > ... > b{2) > b[^) ^ 0, * indicates synrd30l-wise multiplication. + 
indicates symbol-wise addition and mod indicates symbol-wise modular reduction. 

2. The method of claim 1 , wherein k^j. 

3. The method of claim 2, wherein the derivation of constituent^ involves symbol-wise reduction modulo ff. 

4. The method of any one of the preceding claims, wherein p s 2. 

5. The method of any one of the preceding claims, wherein at least two of said constituents identify elements of two 
different respective predetermined sutisets. 

6. The method of any one of the preceding claims, wherein said constituents are derived by applying k respective 
decoding steps. 

• 

7. The method of claim 6, wherein each decoding step yields a set of row coefficients of a generator matrix of a linear 
p-ary code and the constituent identified by that step is derived from the respective set of row coefficients. 

8. The method of claim 6 or daim 7, wherein constituent, is derived by applying a decoding step in respect of said 
input word, constituent is derived by applying a decoding step in respect of said input word modified as a function 
of constituents , and each remaining constituent is derived by applying a respective decoding step in respect of a 
modified input word used for derivation of the preceding constituent, further modified as a function of tliat preceding 
constituent. 

9. The method of any one of claims 6 to 8. wherein each decoding step conprises a single decoding operation. 

1 0. The method of any one of claims 6 to 8, wherein at least one predetermined decoding step comprises a plurality of 
decoding operations. 

1 1 . The method of claim 10, wherein said plurality of decoding operations comprising a decoding step are terminated 
if a parameter evaluated during the decoding operations exceeds a predetermined threshold. 

1 2. The method of daim 1 1 , wherein said threshold is three-quarters of a maximum value attainable by said parameter. 

1 3. The method of any one of daims 9 to 1 2, wherein each said decoding operation involves a Hadamard transform or 
equivalent thereof. 

14. The method of any one of the preceding claims, wherein each symbol of said input word has an integer value. 

15. A method of identifying a code word within a predetermined set of code words, the code words containing symlxsts 
each having 2^ different possible values where j is an integer greater than 1 , comprising the steps of 

receiving an input word; 

deriving the result of redudng the input word modulo 2** ; 

deriving a first constituent, constituenti . by applying a Hadamard transform operation, or equivalent thereof, in 

respect of the result of said reduction modulo 2^ ; 

setting a modifier equal to constituent^ ; 

for each successive value of a counter x from 2 to /: 

deriving the result of symbol-wise sutitraction of tiie current value of the modifier from the input word mod- 
ulo 2'; 

deriving another constituent, constrtuent;^. by applying a Hadamard transform operation, or equivalent 
thereof, in respect of the result of said subtraction modulo 2^; 

setting the modifier equal to the symbol-wise sum modulo 2^ of tiie modifier and tiie product of 2^*'' times 
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constituentj^; 

and 

identifying the code word as having the final value of said modifier. 



16. Apparatus for identifying a code word within a predetermined set of code words, the code words containing sym- 
bols each having pi different possible values where p is a prime and j is an Integer greater than 1 , comprising: 

a decoder for deriving in respect of an input word a plurality of constituents, constituent^ to consfituentj^. of said 
10 code word, each constituent identifying an element of a predetermined subset of a predetermined linear code: 

and 

a combiner for combining said k constituents to identify said code word according to a relationship 



(p^^'^^*constituent^ + p^^'^"^^*constituent;f.^ + ... + p*^^^*constituent2) + p'^^^^*constituenti) mod p^ 

(p*W*constituent;f+p^^'^"^J*constituent/f.i+ ... +p^^^^*constituent2+p^^^^*constituenti) mod pi where each b{x) 
is an integer. J > b{k) > £)(/c-1) > ... > ^(2) > b{^) ^ 0, * indicates symbol-wise muHipiication, + indicates symbol- 
wise addition and mod indicates symbol-wise modular reduction. 
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